Reference: Bio! Hull. 191 : 33-41. (August, 19%) 


Subunit Composition of the Crustacean Hemocyanins: 
Divergence in Incipient Speciation 


CHARLOTTE P. MANGUM AND AMANDA L. McKENNEY* 
Department of Biology, College of William and Mary, Williamsburg , ) Irginia 231S7-8795 


Abstract. The monomeric subunit composition of the 
oxygen carrier hemocyanin was examined in samples of 
the Sesarma retieulatum complex and of Uea minax , 
both of which are believed to be in the process of specia¬ 
tion. The samples were taken on the Atlantic and Gulf of 
Mexico coasts, from disjunct populations that are be¬ 
lieved to have been isolated by the Florida peninsula 
since the retreat of the last glacier. In Atlantic and Gulf 
samples of the S. retieulatum complex, which is believed 
to be in the terminal stages of speciation, the hemocya¬ 
nins differed qualitatively. Several electrophoretic bands 
found in one group, including an invariant band, were 
totally absent from the other. This difference exceeds 
that reported in a previous investigation of a variety of 
polymorphic allozymes in this species complex. It also 
exceeds the physiologically labile differences in hemocy¬ 
anins found previously within a panmictic species of 
brachyuran crustaceans. In U. minax, which is believed 
to have diverged less, Atlantic and Gulf animals ex¬ 
pressed the same number of electrophoretic bands at ex¬ 
actly the same positions. Nonetheless, highly significant 
differences in band frequencies distinguished both Atlan¬ 
tic samples from the Gulf sample, and somewhat less sig¬ 
nificant differences distinguished the Atlantic cold tem¬ 
perate zone samples from the warm temperate zone 
ones. The phenotypes of the major bands, defined as 
those present in high densities, qualitatively distin¬ 
guished Atlantic from Gulf animals, but they did not 
differentiate the two Atlantic samples. The difference be¬ 
tween Atlantic and Gulf members of this species also ex¬ 
ceeds that found previously among polymorphic allo¬ 
zymes. These findings further support the hypothesis 
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that the hemocyanins are among the first proteins to di¬ 
verge structurally in brachyuran speciation. 

Introduction 

With few exceptions, the crustacean hemocyanins 
(Hcs) are composed of many different polypeptide 
chains. When separated by charge in native polyacryl¬ 
amide electrophoresis (PAGE), a dozen or more Cu-con- 
taining bands have been reported in some species (Man- 
gum and Greaves, 1995; Mangum, 1996). In most spe¬ 
cies investigated thus far, the unusual heterogeneity 
provides the basis for intraspecific polymorphisms of a 
magnitude that is unprecedented in the literature on al- 
lozyme variation (Mangum, 1990, 1993a, 1996; Calli- 
cott and Mangum, 1993; Mangum and Greaves, 1996). 
In spite of these polymorphisms, the PAGE banding pat¬ 
tern appears to be species specific, even in sibling and 
previously cryptic species (Reese and Mangum, 1994; 
Mangum, 1996). The specificity suggests that the oxygen 
carrier diverges early in speciation, a hypothesis that is 
consistent with the quantitative importance of crusta¬ 
cean Hcs in aerobic metabolism (Truchot, 1992). 

If the hypothesis is correct, one might expect to find 
differences between populations that are in the process 
of speciation. Below the species level, however, the avail¬ 
able information on structural divergence of the crusta¬ 
cean Hcs appears to be contradictory. On the one hand, 
hybrid lobsters produced in the laboratory by unnatural 
matings of sibling, allopatric parental species can be 
readily distinguished from both parents (Mangum, 
1993b). Less convincingly, because the comparison is 
based on different procedures used in different labora¬ 
tories, samples from two subspecies of a fiddler crab 
differed so much that no bands were clearly shared by the 
two (Mangum and Greaves, 1996). On the other hand, 
no difference was found between the stone crab Menippe 
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mercenaria and the hybrid progeny of postulated natural 
matings between it and its sibling and former cryptor M. 
adina (Mangum, 1996). Also less convincing for the 
same reason described above, the phenotypes found in 
another pair of fiddler crab subspecies were quite similar 
(Mangum and Greaves, 1995). 

In the present investigation, we analyzed the banding 
patterns of He monomers found in populations isolated 
by the Florida peninsula since the retreat of the last gla¬ 
cier. One pair is believed to be in the terminal stages of 
speciation, and the other group is believed to have di¬ 
verged less. Because each sample proved to be polymor¬ 
phic, we also characterized the intraspecific variation in 
detail, to ensure against confusing it with geographic di¬ 
vergence. 

The species investigated are (1) the grapsid marsh crab 
Sesarma reticulatum (Say) from the Atlantic coast of 
North America and its sibling from the Gulf of Mexico, 
which has not yet been described and is currently known 
as Sesarma sp. (near reticulation) (Zimmerman and 
Felder, 1991), and (2) the fiddler crab Uca mimix 
(LeConte), also from both coasts. 

Materials and Methods 


Animals 
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I igure 1 . 1 he frequency distribution of the total number of bands 

found in an individual of Sesarma reticulatum (n = 97) and Sesarma 
sp. (// = 64). 


their monomers and thus minimizes light scattering. He 
concentration was estimated from the absorbance of 
monomers, using the extinction coefficient for Cancer 
borealis He (Nickerson and van Holde, 1971), and addi¬ 
tional aliquots of material were then diluted with the 
same butter to produce a concentration of about 
5 og l . We emphasize the importance of the dissocia¬ 
tion step in this procedure, which, in our experience, 
eliminates an error of about 15%. 


All individuals investigated appeared to be intermolt 
crabs, as indicated by hard exoskeletons. 

Atlantic populations of Sesarma reticulatum were 
sampled at Taskinas Creek and Indian Field Creek, both 
of which empty into the York River estuary in Virginia. 
Gulf of Mexico populations of Sesarma sp. were sam¬ 
pled at Cocodrie and Joseph's Harbor (Rockefeller ref¬ 
uge) in Louisiana, and at Ocean Springs in Mississippi. 
Atlantic cold temperate zone populations of Uca mina.x 
were sampled at five sites on the York River estuary 
(Ware Creek, Croaker Landing, York River State Park, 
Kings Creek, and Indian Field Creek) and one on the 
James River estuary (College Creek) in Virginia; these 
sampling sites are distributed along a salinity gradient 
ranging from oligohaline to upper mesohaline. A warm 
temperate zone population was sampled at Holden 
Beach in North Carolina. Gulf populations were sam¬ 
pled at Cypremort Point in Louisiana and Fort Bayou 
near Ocean Springs in Mississippi. 

Preparation of material 

Animals were bled immediately with a hypodermic sy¬ 
ringe, and the blood was frozen until shortly before 
analysis. After thawing, serum was expressed from the 
clot by centrifugation, and an aliquot was diluted with a 
buffer (50 mmol l _l Tris HC1, pH 8.9, containing 
10 mmol 1 1 EDTA) that dissociates He oligomers to 


Electrophoresis 

Native PAGE of dissociated material was performed 
at constant current, using the procedures and reagents 
detailed by Hames and Rickwood (1990). Two volumes 
(10 and 3 /j 1) of the material from each individual were 
examined. Most gels were stained only with Coomassie 
blue. To detect Cu, however, examples from each sample 
that collectively expressed all bands found in the sample 
were first electrophoresed, and then inspected on a long¬ 
wave UV transilluminator in the presence of bathocu- 
proine sulfonate (Bruyninckx et aL, 1979). To estimate 
the proportions of the bands, examples were scanned 
with a BioRad Video Densitometer (Model 300A), and 
the scans were quantified with the ID Analyst software 
package. 

Results 

Sesarma reticulatum complex 

No difference was found between the individuals col¬ 
lected in Louisiana (n = 61) and Mississippi (3), and none 
was found between the individuals collected from the two 
tidal creeks in Virginia (n = 97). In both cases the data 
were combined for analysis and presentation. All bands 
were positive for Cu, and no band was related to sex. 

Whereas all but a few individuals in the Atlantic sam¬ 
ple expressed 6 He bands, the modal number in the Gulf 
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Sesarma reticulatum (Atlantic) 



Sesarma sp. (Gulf) 



Band No. 

Figure 2. The frequencies of each band in Sesarma reticulatum and 
Sesarma sp. Sample size as in Fig. 1. 


sample was 7 (Fig. 1). A total of only 9 bands was found 
in the Atlantic animals, but a total of 1 1 occurred in the 
Gulf animals (Figs. 2 and 3). Bands 2, 4, 5, and 11-13 in 
the two groups co-migrated to identical positions (Fig. 
3). Two of these, however, differed in variability. Band 4 

Sesarma sp. 
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was invariant in Sesarma sp. but not in S. reticulatum 
(Fig. 2), in which it was often found in lower density as 
well (Fig. 3). Band 12 was variable in Sesarma sp. but 
not in S. reticulatum (Fig. 2). 

In both groups a cathodic cluster of bands separated at 
some distance from an anodic cluster (Fig. 3). In both 
groups, the most common phenotype included bands in 
each cluster that co-migrated with bands in the other 
group, namely bands 2, 4, and 11-13. Nonetheless, the 
two groups were reliably distinguished by the presence in 
Sesarma sp. of a pair of both frequent and dense bands 
(8 and 9) in between the two clusters; because it was in¬ 
variant, band 9 was an absolute marker. Members of this 
pair of bands were never found in S. reticulatum (Figs. 2 
and 3). Additional, though less striking, differences in¬ 
clude (1) one or two cathodic bands in each group that 
were not found in the other (1 in Sesarma sp. and 6 in 
S. reticulatum ), and (2) two infrequent and less dense 
members (7 and 10) of the set of bands in between the 
cathodic and anodic clusters in Sesarma sp. Band 14 was 
found only once and only in trace densities in Sesarma 
reticulatum; it is unlikely to be a useful marker. 

Both groups were polymorphic: 7 phenotypes were 
found in S . reticulatum, and 13 were found in Sesarma 
sp. (Table I). Of the 1 1 bands in Sesarma sp., 6 were vari¬ 
able, and the 5 qualitatively invariant bands still varied 
in density (Fig. 3). Although only 4 of the 9 bands were 
qualitatively variable in S . reticulatum (Fig. 2), all but 
one of the invariant bands varied in density (Fig. 3). 

With the exception of the infrequent band 14, the vari¬ 
ation in S . reticulatum was limited to bands 3, 4, and 5, 
which varied concomitantly in the classic pattern ex¬ 
pected of the products of alleles encoded at the same lo- 
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Figure 3. Diagrammatic representation of the electrophoretic bands in Sesarma. If no range is speci¬ 
fied, then band density did not vary. 
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Table I 


Phenotypes observed in Sesarma 


A. Sesarma retieulatum 
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cus. Each of these three bands was expressed either alone 
or together with one of the other two, but all three were 
never found together (Table 1); nor were all three ever 
absent coincidentally. In general, density was highest 
when one of the three occurred alone (Fig. 4). Of the 


Sesarma sp 


A BC 



D E F 



members of the triplet, band 4 was by far the most fre¬ 
quent (Fig. 2). In Sesarma sp., band 4 was invariant, and 
its density did not clearly decrease when band 5 was ex¬ 
pressed. The third member ol the S retieulatum triplet 
(band 3) was never found in Sesarma sp. 

In Sesarma sp., bands 8 and 12 also appeared to vary 
as expected of the products of a single locus. Again, each 
was found alone or together with the other, but one or 
both were always expressed. In this case, the variation in 
density (Fig. 3) appeared to be independent. The fre¬ 
quencies of each were quite similar in Sesarma sp., but 
band 8 was never found in S . retieulatum (Fig. 2). 

Because this is the only of now several investigations 
in which we have encountered allozyme-like variation, 
we have chosen not to estimate "gene" frequencies and 
Hardy-Weinberg ratios until the genetic character of the 
observed variation can be corroborated. We do, how¬ 
ever, report the frequencies of each phenotype in the Ses¬ 
arma complex (Table 1). 

Uca minax 
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Figure 4. Examples of the variation in Sesarma. With the excep¬ 
tion of F, each panel shows higher (bottom) and lower (top) concentra¬ 
tions of the material from the same individual. 


No difference was found between the Atlantic cold 
temperate zone animals (/? = 161) collected at the several 
sites in Virginia; therefore, they are treated here as a sin¬ 
gle sample. Similarly, no difference was found between 
the 22 Gulf coast animals collected in Louisiana and the 
49 collected in Mississippi, which are also treated here as 
a single sample. In each area, that is. Gulf coast and At¬ 
lantic cold (Virginia) temperate zone, a total of 17 bands 
were found. In any one sample, each band co-migrated 
to the same position as did a counterpart in the other two 
samples (Fig. 5). In other words, no marker bands were 
diagnostic of a particular sample. In addition, the modal 
number of bands expressed in each sample is the same 
(Fig. 6). Nonetheless, each sample differs significantly 
from the other two. 
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Figure 5. Diagrammatic representation of the electrophoretic 
bands in Vca minax If no range is specified, then band density did not 
vary'. 


First, the frequencies of the bands differ (Fig. 7). Pair¬ 
wise comparisons were made in a contingency test of in¬ 
dependence. Because some cells in the Atlantic warm 
temperate zone sample were empty, the Chi-square sta¬ 
tistic was used. Each of the two Atlantic samples differs 
from the Gulf sample at a probability level several orders 
of magnitude less than P = 0.001. Second, the variation 
in band density was similar in the two Atlantic samples, 
but it differed from that in the Gulf sample (Fig. 5). Per¬ 
haps most conspicuous, the phenotypes of the major 
bands, defined as those that comprise the majority of the 
material present, qualitatively distinguished the two At¬ 
lantic samples from the Gulf sample (Table 11). Although 
the phenotypes of the major bands in the two Atlantic 
samples were the same, their frequencies differ signifi¬ 
cantly. Again, in a pairwise comparison by a contingency 
test of independence, the two differ at about P = 0.001. 

Each of the three samples exhibited considerable poly¬ 


morphism, with a very large number of phenotypes (Ta¬ 
ble 111, Fig. 8). The sample from the Gulf coast, the small¬ 
est of the three (/? = 71), had the smallest number of phe¬ 
notypes (26). The sample from the warm temperate zone 
on the Atlantic coast, only slightly larger (n = 73), had a 
similar number of phenotypes (28). The sample from the 
cold temperate latitudes on the Atlantic coast, which was 
the largest (/? = 161), exhibited a proportionately larger 
number of phenotypes (56), which is somewhat surpris¬ 
ing. One would expect the rate of increase to diminish 
with increasing sample size. 

In all three samples, the bands that occurred in either 
high densities or high frequencies were clearly positive 
for Cu. Bands 1,2, 6, and 7 were neither clearly positive 
nor clearly negative because the combination of low den¬ 
sity and low frequency precluded repeatedly overloading 
the gels for the less sensitive Cu detection procedure. 

The only band that was invariant in all three samples 
was 14 (Fig. 8). However, bands 3, 4, and 13 were absent 
in only a very few individuals, and the variation in band 
16 was confined to the Gulf sample (Fig. 7, Table 111). 

None of the qualitative variation of any one band was 
coupled to that of any other in a way that clearly suggests 
encoding at a common locus (Table 111). One instance 
of coupled quantitative variation was noted: in the two 
Atlantic samples, the densities of bands 9 and 11 always 
varied concomitantly and directly. 

Discussion 

On the basis of allozymic variation at 13 polymorphic 
loci, Felder and Staton (1994) concluded that the genetic 
distance between Atlantic and Gulf of Mexico samples 
of Sesartnci reticulatum (Say) is similar to that found in 
recently speciated populations of other brachyurans. 
The data strongly supported the authors' earlier conclu¬ 
sions, which were based on coloration, osmotic re¬ 
sponses, and reproductive chronology (Zimmerman and 
Felder, 1991: Staton and Felder, 1992). All allozyme 
differences between samples of the two disjunct popula- 
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Figure 6. The frequency distribution of the total numher of bands found in an individual of Vca minax. 
Values of n are 161 for the cold temperate zone sample, 73 for the warm temperate sample, and 7 1 for the 
Gulf sample. 
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Atlantic 

Cold temperate (solid) vs. warm temperate (hatched) 



Gulf (solid) vs. Atlantic warm temperate (hatched) 



Figure 7. The frequencies of each band in i ca minax. Sample sizes 
as in Fig. 5. 

lions, which they designated S. rcliatlantm and Sesamia 
sp., were quantitative, although several morphs com¬ 
pletely absent from one coast existed in very high fre¬ 
quencies on the other. 

The Hcs of the two Scsarma samples investigated here 
dilfer qualitatively. Band 9, found in all individuals of 
Scsarma sp.. was totally absent from the N. retiathUum 
sample. Therefore, the present findings not only support 
the hypothesis of speciation, they also suggest that the 
Hcs differ more than the 13 polymorphic enzymes 
(much less an additional 7 enzymes, which were mono- 
morphic in both species). This suggestion is constrained, 
however, by the use of different sampling sites on the At¬ 
lantic coast in the two investigations. Whereas Felder 
and Staton (1994) examined only warm temperate zone 
members of S’, relicalatwn, we examined only cold tem¬ 
perate zone representatives. Populations in the two bio¬ 
geographic provinces on the Atlantic coast are by no 
means disjunct, but Cape Hatteras, which separates 
them, is a major faunal discontinuity (c.g., Vernbergand 
Vernberg, 1972) and, in several taxa, populations found 
on either side of it differ (c.g.. Vca minax in the present 
investigation). 

Felder and Staton (1994) reached different conclu¬ 
sions for Uca minax. In this ease the allozyme variation 
was deemed insufficient to warrant specific separation. 
Although the inference was also based on 13 polymor¬ 
phic loci, the variation at 9 of them was quite small. The 


only pertinent morphological evidence in this species is 
the intensity of the red coloration on the articles of the 
ehelipeds and legs, which is strongly developed in Atlan¬ 
tic and western Gulf populations but, oddly, not in east¬ 
ern Gulf populations. No comparative information on 
physiological or reproductive responses is available. As 
indicated above, our samples of eastern and western Gulf 
members of V minax did not differ, indicating no rela¬ 
tionship to the color patterns. Most of our sample origi¬ 
nated from a Mississippi site where the color difference 
from Atlantic animals is maximal, and the remainder 
was collected in Louisiana, where the pattern more 
closely resembles that found in Atlantic members of the 
species. A similar conclusion was reached by Felder and 
Staton (1994), v\hose dendrogram based on allozymic 
differences does not reflect the color patterns. 

Regardless, we still found clear, though primarily 
quantitative, differences in the overall banding patterns 
of the He monomers in each of the three samples ana¬ 
lyzed. The same number of Cu-containing bands was 
found in Atlantic and Gulf animals, and each band in 
one sample co-migrated with a counterpart in the other. 
The most probable inference is that the co-migrants rep¬ 
resent the same polv peptides. Nonetheless, the frequen¬ 
cies of the bands are very different in each sample. The 
overall differences between the two Atlantic samples are 
no greater than those observed within Callincctes sapi- 
chts , which are not genetically fixed in the adult stage 
(Mangum, 1990; deFur cl al., 1990). In contrast, the 
differences between either Atlantic sample and the Gull 
sample are clearly greater than those in C. sapidus. which 
can be induced by environmental change. Although 
cross-acclimation experiments should be performed in 


Table 11 

Frequencies of I fie phenotypes of the major hands . defined as those that 
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Tabic 111 


Phenotypes observed in Uca minax 
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FigurcS. Examples of lhe variation in Vea mtnax Top: Atlantic animals. With the exception of panel 
F, lhe lanes are arranged in pairs. Each pair of lanes, from left to right, shows a higher and a lower concen¬ 
tration of material from lhe same individual. Bottom: Gulf of Mexico animals. All lanes are arranged in 
pairs showing two concentrations of material from the same individual. With the exception of lhe pair at 
the far left, which shows first the lower and then lhe higher concentration, each pair alternates first lhe 
higher and then the lower concentration, as above. 


V . minax, we suggest tentatively that the differences may 
reflect, at least in part, genetic divergence of the disjunct 
Atlantic and Gulf populations. This suggestion, if true, 
also supports the hypothesis that the Hcs diverge early in 
speciation. 

We note that the differing magnitudes of the diver¬ 
gence in the two complexes examined here agree with the 
findings of Avise and co-workers as well as with those of 
Felder and students. On the basis of results for 19 species 
ranging from bivalves to chelicerates to teleosts, Avise 
(1992) found no evidence of a uniform molecular clock 
that measures the rate of divergence of Gulf and Atlantic 
populations. 

We also suggest tentatively that very early divergence 
of the Hcs may be accelerated when the products of spe¬ 
ciation must adapt to different thermal regimes (see Re¬ 
ese and Mangum, 1994), such as those found in the 


warm temperate (Gulf) and cold temperate (Atlantic, 
north of Cape Hatteras) zones. 

Perhaps most important from a physiological point of 
view, the Atlantic and Gulf samples of U. minax are 
qualitatively differentiated by the phenotypes of the ma¬ 
jor bands, those that are most likely to influence respira¬ 
tory properties. In addition, this inference remains valid 
when the comparison is restricted to Atlantic warm tem¬ 
perate zone vs. Gulf samples. The thermal regimes expe¬ 
rienced by these two populations differ relatively little, 
certainly far less than the thermal regime experienced by 
either population and that experienced by Atlantic cold 
temperate zone populations. 

Deevev (1950) suggested that environmental tempera¬ 
ture is the primary factor responsible for the disjunction 
of Atlantic and Gulf species, a widespread phenomenon 
in a number of animal phyla (see also Avise, 1992). In 
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his scenario, disjunct populations are relicts of the retreat 
of the last glacier. The concomitant rise in temperature 
on both coasts around the Florida peninsula caused the 
retreat northwards of temperate zone species and their 
replacement in Florida by tropical and subtropical spe¬ 
cies invading from the south. Barnwell and Thurman 
(1984) correctly pointed out that biological factors such 
as the attendant changes in vegetation must also play a 
role in the disjunction of fiddler crab species. Whereas U. 
minax, for example, inhabits Spartina salt marshes, the 
species of Uca now found in subtropical regions of the 
Florida peninsula are mangrove dwellers. 

Although we do not doubt the importance of habitat, 
we believe that temperature is likely to be extremely im¬ 
portant in the present case. Temperature is the factor 
most often associated with divergence of respiratory 
properties among the Hcs of closely related species 
(Reese and Mangum, 1994). Consequently, the thermal 
sensitivities of the respiratory properties of Hcs express¬ 
ing the most frequent phenotypes of major bands in Uca 
niina vand the two Sesanna species will be of consider¬ 
able interest. 
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